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Abstract 

We present a detailed study of direct CP violation and branching ratios in the 
channels B 0,ziz — ► 7r + ir~V 0,zil , where V is a vector meson (K* 0,:iz or p^). Emphasis 
is placed upon the important role played by p° — lo mixing effects in the estimation 
of the CP-violating asymmetry parameter, a cp , associated with the difference of B 
and B decay amplitudes. A thorough study of the helicity amplitudes is presented 
as a function of the pion-pion invariant mass. All of the calculations and simulations 
considered correspond to channels which will be analyzed at the LHCb facility. 
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1 Introduction 



Understanding the physical origin of the violation of CP (Charge Conjugation x Parity) 
symmetry is one of the main goals of Particle Physics at the present time. Recent ex- 
periments at e + e~ colliders (BaBar, Belle) have produced fundamental results which 
strengthen the CKM picture of CP violation QjJ |2] in the B meson sector [HUE]. How- 
ever, the main results of these two collaborations are related to B decays into pairs of 
pseudo-scalar mesons or into a vector plus a pseudo-scalar meson. 

A very broad physics program can also be carried out in the sector with two vector 
mesons in the final state, following B decay. Apart from measuring the standard angles, 
a, (3 and 7 of the Unitary Triangle (UT), the vector mesons are polarized and their decay 
products (usually long-lived CU + mesons) are correlated. This opens the possibility of 
making interesting cross-checks of the Standard Model predictions as well testing some 
specific models of CP violation beyond the SM approach (BSM). 

In the special case of two neutral vector mesons, the orbital angular momentum, £, 
and the total spin, S, of the V®V 2 ° system satisfy the equality t = S = 0, 1,2. The CP 
eigenvalues are defined as (— l) e . Because of the allowed values for the angular momentum, 
£, one has a very clear indication of any mixing of different CP eigenstates and hence of 
CP non-conservation. The separation of the different CP eigenstates requires a detailed 
analysis of the final angular distributions [5 . However, because this analysis can be 
carried out in a model-independent way, it provides a significant constraint on any model. 

After explaining the helicity formalism (Section II), a special study is devoted to the 
final state interactions (FSI) and the key role of p° — uo mixing (Section III). A complete 
and realistic determination of the helicity amplitudes, in the framework of the effective 
Hamiltonian approach, is introduced in Section IV. Then, the main results of the Monte- 
Carlo simulations, providing estimates of the various density matrix elements hij, are 
shown in Section V. In the following section (Section VI) the numerical analysis and 
discussions about the branching ratios and asymmetries for B decays into two vectors 
(B -> p°(uj)V 2 , with V 2 = K*°,K*°,K*-,K* + ,p + ,p-) are given in detail. These two 
vectors, p°(cu) and V 2 , each decay into two pseudo-scalars. Emphasis is put on the angular 
distributions of the pseudo-scalar mesons in both the helicity and transversity frames. 
Finally, in the last section, we summarize our results for the different channels which will 
be investigated in future experiments at pp colliders and make some concluding remarks. 

2 General formalism for B — > V1V2 decays 
2.1 Helicity frame 

Because the B meson has spin 0, the final two vector mesons, Vy and V 2 , have the same 
helicity, Ai = A2 = —1,0,-1-1, and their angular distribution is isotropic in the B rest 
frame. Let H w be the weak Hamiltonian which governs the B decays. Any transition 
amplitude between the initial and final states will have the following form: 

H X =(V 1 (X)V 2 (X)\H W \B) , (1) 

where the common helicity is A = —1, 0, +1. Then, each vector meson Vi will decay into 
two pseudo-scalar mesons, a*, b{, where Oj and foj can be either a pion or a kaon, and the 
angular distributions of a, and 6, depend on the polarization of Vi. 



The helicity frame of a vector-meson Vi is defined in the B rest frame such that the 
direction of the Z-axis is given by its momentum, fy. Schematically, the whole process 
gets the form: 

B — > Vi + V 2 — ► (at + 61) + (a 2 + b 2 ) . 

The corresponding decay amplitude, M X (B — > Yli=i( a i + i s factorized out according 
to the relation, 

2 2 

M a( 5 ~> Z>* + = ^a(5 ^ + V 2 ) x JJ -> fli + k) , (2) 

i=i i=i 

where the amplitudes Ai(Vi — > a>i + bj) are related to the decay of the resonances Vi. The 
Ai(Vi — > a,i + 6j) are given by the following expressions: 

1 

mi=— 1 
1 

A 2 (F 2 ^a 2 + fc 2 ) = ]T c 2 Dl m2 (<p,8 2 ,-<P) . (3) 

m 2 =— 1 

These equalities are an illustration of the Wigner-Eckart theorem. In Eq. (|3*|). the ci and 
c 2 parameters represent, respectively, the dynamical decays of the Vi and V 2 resonances. 
The term D\ m .((j)i,6i, —(pi) is the Wigner rotation matrix element for a spin-1 particle 
and we let A(aj) and A(6j) be the respective helicities of the final particles a, and fej in 
the Vi rest frame. #1 is the polar angle of a\ in the V\ helicity frame. The decay plane of 
V\ is identified with the (X-Z) plane and consequently the azimuthal angle <px is set to 0. 
Similarly, 6 2 and <fi are respectively the polar and azimuthal angles of particle a 2 in the 
V 2 helicity frame. Finally, the coefficients rrii are defined as: 

rrii = A(0i) - A(6i) . (4) 

Our convention for the D\ m .(a, f3, 7) matrix element is given in Rose's book [6 , namely: 

D x, mi ( a ^n) = exp[-i(Aa + ma)} 4,™^) • ( 5 ) 

The most general form of the decay amplitude M. [B — > Yli=i( a i + M) * s a linear super- 
position of the previous amplitudes M\(B — > X^=i( a « + M) denoted by, 

2 2 
^( S _>^( a . + 6 .)) =^M A (fi^^(a i + 6 i )) • (6) 



i=l A i=l 



The decay width, T(5 — > V1V2), can be computed by taking the square of the modulus, 



M(B ^ ElM + h)) 



, which involves the three kinematic parameters 8i,6 2 and 4>. 



This leads to the following general expression: 

1 2 2 

A i=l A, A' 



which involves three density-matrices, h\ t y, F\ t y(9i) and G\ t y(9 2 , 0)- The factor h\ t y = 
H\Hy is an element of the density-matrix related to the B decay, while F\ t \>(Qi) represents 
the density-matrix of the decay V\ — > a\ + b\. In a similar way, G\^{9 2 , 0) represents the 
density- matrix of the decay V 2 — > a 2 + b 2 . 

The analytic expression in Eq. (|7J) exhibits a very general form. It depends on neither 
the specific nature of the intermediate resonances nor their decay modes (except for the 
spin of the final particles). This approach also presents three key advantages. The first 
one comes from the fact that all the dynamics of the B decay is introduced into the 
coefficients h\ t \i. This allows us to use various theoretical models involving different 
dynamical processes and form factors. The second is associated with the fact that the 
formal expressions for F\ t \>(9i) and G\^y{9 2 , 0), which are related to the polarization of 
the intermediate resonances, remain unchanged whatever the coefficients /ia,a' happen to 
be. Finally, correlations among final particles arise in a straightforward way because of 
the previous expression which relates the angles 9i, 9 2 and 0. Consequently, a probability 
density function (pdf) can be inferred from the general expression and one gets: 

f(9 9 rh) - ^ - VlV2) (R) 

A 2 ' 9) T(B -> ViV 2 )d(cos9 1 )d(cose 2 )d<l> ' { } 

where the angles 9\,9 2 and were defined earlier and T(B — > ViV 2 ) is the partial decay 
width. This function allows one to compute three other pdfs separately for the variables 
0i, 02 and 0. 

The previous calculations are illustrated by the reaction B° — ► K*°p° where K*° — > 
K + 7i~ and p° — > 7r + 7r _ . In this channel, since all the final particles have spin zero, the 
coefficients mi and m-i, defined in Eq. (|3J), are equal to zero. The three- fold differential 
width has the following form: 



d 3 T{B ^V^) 



oc (h ++ + h )sin 2 6 ) isin 2 6 ) 2/4 + /i ocos 2 6 , 1 cos 2 6 ) 2 



c?(cos 6 l 1 )<i(cos 02)d(p 

+ j^fe(/i +0 ) cos 4> — ,ym(h +0 )sm <p + &e{h J)cos <p — J%i(/i _)sin ^jsin 26 l 1 sin 2^2/4 

+ ji^O + „)cos 20 - J%i(/i + _)sin 20|sin 2 ^isin 2 ^ 2 /2 , (9) 

where all the terms in Eq. Q have been already specified. It is worth noticing that the 
expression in Eq. (JHJ) is completely symmetric in 0\ and #2 and consequently, the angular 
distributions of a\ in the V\ frame is identical that of a 2 in the V 2 frame. From Eqs. (jSji 
and © the normalized pdfs of 0\, 2 and can be derived and one finds: 

/(cos0 lj2 ) = (3h 00 - l)cos 2 6»i >2 + (1 - /ioo) , 

g(<t>) = 1 + 2 ^(/i + _)cos 20-2 J^m(/i + _)sin 20 . (10) 



2.2 Transversity frame 

Initially, the transversity frame (TF) was introduced by A. Bohr [7j in order to facilitate 
the determination of the spin and parity of a resonance decaying into stable particles. It 
can be extended to a system of two vector mesons coming from a heavy meson, B or B, in 
order to perform tests of CP symmetry. In displaying new angular distributions, the TF 



provides complementary physical information to that seen in the standard helicity frame. 
The construction of the TF and its use require several steps. For a clear illustration, see 
Fig. H where the channel B° — > p°K*° is chosen. 

Departing from the B rest frame, the common helicity axis, (Ah), is given by the 
direction of the momentum p[. This and the decay plane, (ITd), of the vector meson 
(K*° — > K + 7r~) are the main ingredients of the TF. The vector meson p° is taken at rest 
(origin of the frame) and the X-axis is given by (A#). In the decay plane, (rip), the 
Y-axis which is orthogonal to the X-axis, is chosen in such a way that the K + meson has 
the Y-component of its momentum greater than or equal to zero. The Z-axis, which is 
orthogonal to the plane {Hd), is obtained by the classical relation ez = eyXey. 

The angular distributions of the ir^ coming from the p° decay are referred to the new 
Z-axis. It is worthy noticing that, in the TF, the flying meson and its decay products are 
very energetic compared to the B frame. Explicitly, the p° energy is given by the relation, 

E p0 = (m B 2 - m x 2 - m 2 2 )/2m l « 17 GeV , (11) 

where m,\ and m,2 refer to the masses of the K*° and p° resonances, respectively. As far as 
the transition amplitudes in the TF are concerned, they are a simple linear combination 
of the helicity amplitudes, namely: 

Hp = ' a,ld Ht = !k 7T L • (12) 

while Hq remains unchanged. We can rewrite the angular distributions given in Eq. |T0|) 
by using the relations from Eq. (fT2|) and angles 6^2, <fi expressed in the transversity frame. 
Thus one gets, 

/ r (cos0i, 2 ) = (3|# T | 2 - l)cos 2 # li2 + (l - \H T \ 2 ) , 

gr(<f>)= (l + \H \ 2 -\H P \ 2 )cos2<p . (13) 



3 Final state interactions and p — uj mixing 

3.1 Factorization hypothesis 

Final state interactions (FSI) represent unavoidable effects in hadronic physics and they 
play a crucial role in heavy resonance decays jH] . In the case of a B meson, characterized 
by a center-of-mass energy ^Js ~ 5.3 GeV, the charmless weak decays of the b-quark lead 
to light energetic quarks which can exchange several gluons amongst themselves as well 
as with the spectator quark in the B meson. This fundamental process occurs in decays 
described by tree, penguin and annihilation diagrams and is characterized by two regimes: 
perturbative and non-perturbative. In order to handle the FSI in both regimes, the usual 
method is inspired by the effective Lagrangian approach. Perturbative calculations at 
next-to-leading order (NLO) are performed for a scale higher than mj (since our analysis 
is focused on B decays) and the non-perturbative effects are inserted for a scale lower 
than mb- This general method is called the factorization procedure [H] and further details 
are given below. 

In the factorization approximation, either the vector meson p°(ui) or the K* is gener- 
ated by one current in the effective Hamiltonian which has the appropriate quantum 



numbers. For the B decay processes considered here, two kinds of matrix element 
products are involved after factorization (i.e. omitting Dirac matrices and color labels): 
(p°(u)\(uu)\0)(K*\(q i q j )\B ± '°) and {K%q iqj )\$) (p Q (ui)\{ub)\B ±fi ) , where q { and y could 
be either u, s or d. We will calculate them in two phenomenological quark models. 

The matrix elements for B — > X* (where X* denotes a vector meson) can be decom- 
posed as follows [TU] . 



(X*\J^\B) = ■ e^ up(T e w f B p^V{k 2 ) + il e*Jm B + m x *)A 1 (k 2 ) 

m B + m x * I F 

"—^{Pb + Px*)»A 2 {k 2 ) - 6 -^2m x , ■ k^A 3 (k 2 " 

tub + mx* k 

+ i e -^2m x * ■ k^A (k 2 ) , (14) 

where is the weak current, defined as J M = 97^(1 — 75)6 with q = u,d,s and 
k = pb ~Px* an d e M is the polarization vector of X*. The form factors A , Ai, A 2 , A 3 and 
V describe the transition 0" — > 1 _ . Finally, in order to cancel the poles at q 2 = 0, the 
form factors respect the condition: 

A 3 (0)=A (0), (15) 

and they also satisfy the following relations: 

A />2x m B + m x * 2 m B -m x * , 2 

A 3 (A; ) = } " 2m x . ^ } " (16) 

In the evaluation of matrix elements, the effective number of colors, N%" , enters through 
a Fierz transformation. In general, for an operator Oj, one can write, 

1 1 +& ,withi = I,--- ,10 , (17) 



(iV c e// ), 3 

here £3 describes non-factorizable effects. £j is assumed to be universal for all the operators 
Oj. Naive factorization assumes that we can replace -in a heavy quark decay- the matrix 
element of a four fermion operator by the product of the matrix elements of two currents. 
This reduces to the product of a form factor and a decay constant. This assumption is 
only rigorously justified at large values of iV c . But it is known that naive factorization 
may give a good estimate of the magnitude of the B decay amplitude in many cases jTT] . 



3.2 FSI at the quark level: strong phase generated by the pen- 
guin diagrams 

Let A be the amplitude for the decay B — > p°(u)K* — > tt + tt^K* (a similar procedure 
applies in the case where we have a p ^2] instead of the K*), then one has, 

A = (K*tt-tt + \H t \B) + (K*tt-tt + \H p \B) , (18) 



with H T and H p being the Hamiltonians for the tree and penguin operators. We can 
define the relative magnitude and phases between these two contributions as follows, 



A = (K*tt-tt + \H t \B)[1 + re i5 e^] , 

A = (K*7i + n-\H T \B)[l + re iS e-^} , (19) 

where 5 and <fi are strong and weak phases, respectively. The phase arises from the 
appropriate combination of CKM matrix elements, = a.ig[(V t bV t * s ) / (V u bV* s )}. As a result, 
sin0 is equal to sin 7, with 7 defined in the standard way The parameter, r, is the 
absolute value of the ratio of tree and penguin amplitudes: 



r = 



{p°(u)K*\H p \B) 



{p°(uj)K*\H T \B} 



(20) 



3.3 Strong phase generated by the p° — oj mixing 

In the vector meson dominance model pi] , the photon propagator is dressed by coupling 
to vector mesons. From this, the p° — uj mixing mechanism ^H] was developed. In order to 
obtain a large signal for direct CP violation, we need some mechanism to make both sin 5 
and r large. We stress that p° — uj mixing has the dual advantages that the strong phase 
difference is large (passing through 90° at the uj resonance) and well known [T2*l ll6|. With 
this mechanism, to first order in isospin violation, we have the following results when the 
invariant mass of the 7t + tt~ pair is near the mass of the u resonance, 

(K*n-n + \H T \B) = -^-U^ + ^t p , 

$p$tjj $p 

(K*n-n + \H p \B) = -^-ft^ + % p . (21) 

Here t v (V = p or uj) is the tree amplitude and py the penguin amplitude for producing 
a vector meson, V, g p is the coupling for p° — ► 7r + 7r~, Yl puJ is the effective p° — uj mixing 
amplitude, and sy is the inverse propagator of the vector meson V, 

sy = s — my + imyYy , (22) 

with yfs being the invariant mass of the 7r + 7r~ pair. We stress that the direct coupling 
uj — > 7r + 7r~ is effectively absorbed into fi^ [T7j, leading to the explicit s dependence of 
Ii pu) . Making the expansion Ii piJJ (s) = n pw (m^) + (s — m^)n^(mf,), the p° — uj mixing 
parameters were determined in the fit of Gardner and O'Connell [TH]: S^e 11^ (m^) = 
-3500 ± 300 MeV 2 , Jm U pu) (m 2 J = -300 ± 300 MeV 2 , and fi^(m 2 ) = 0.03 ± 0.04. In 
practice, the effect of the derivative term is negligible. From Eqs. (fTH|) and (|21jl one has 

re^e^ = ^ pujPuJ + SujPp (23) 



Defining 



— = r'e l ^ +<t>) , - = ae iS " , ^ = (3e i& ? , (24) 
tp tp 



where 5 a , 8p and S q are partial strong phases (absorptive part) arising from the tree and 
penguin diagram contributions. Substituting Eq. (J24"j) into Eq. (|23|) . one finds: 



re = re '— — J — , (25) 

where the total strong phase, 5, is mainly proportional to the ratio of the penguin and 
tree diagram contributions. 



3.4 Importance of the strong phase for BB asymmetry 

Under a CP transformation the strong phase, 5, remains unchanged, while the weak phase, 
4>, which is related to the CKM matrix elements, changes sign. Thus, the asymmetry 
parameter, a^p, which can reveal direct CP violation, can be deduced in the following 
way: 

dir = A 2 - A 2 = -2 r sm5 sin0 

CP A 2 + A 2 l + r 2 + 2r cos 5 cos ' 1 1 

It is straightforward to see that the parameter a^ p depends on both the strong phase 
and the weak phase and, consequently, that the maximum value of a^p can be reached 
if sin 5 = 1. This is why the strong final state interaction (FSI) among pions coming 
from p° — uj mixing enhances the direct CP violation in the vicinity of the mass of the uj 
resonance. 

In the Wolfenstein parametrization jTH], the weak phase comes from [V^V^/K^KTJ 
and one has for the decay B — > p (u)K*, 

-v 



sm . 



Vp 2 + v 



2 



cos0 = = = , (27) 
V P 2 + V 2 

while the weak phase comes from [VtbV t * d /V u bV* d ] for the decay B — > p°(u))p, 

r, 



sm ■ 



y/\p{\- p) -r, 2 ] 2 + ri 2 



cos0= - ^-^ 2 , ( 2 f 

v / [p( l_ p) _ r? 2 ]2+ , / 2 



The values used for p and r\ will be discussed in Section V. 



4 Explicit calculations according to the effective 
Hamiltonian 

4.1 Generalities concerning the OPE for weak hadronic decays 
4.1.1 Operator product expansion 

The operator product expansion (OPE) [20] is an extremely useful tool in the analysis of 
weak interaction processes involving quarks. Defining the decay amplitude AiM —>■ F) 



as 



A{M^F) <xCi(fi)(F\Oi(n)\M) , 



(29) 



where Cj(^) are the Wilson coefficients (see Section R. 1.2)1 . Oj(/x) are the operators given 
by the OPE and \i is an energy scale, one sees that the OPE separates the calculation of the 
amplitude, AiM — > F), into two distinct physical regimes. One is related to hard or short- 
distance physics, represented by Cj(/z) and calculated by a perturbative approach. The 
other is the soft or long-distance regime. This part must be treated by non-perturbative 
approaches such as the 1/N expansion [21], QCD sum rules |22j . hadronic sum rules or 
lattice QCD. 

The operators, Oj, are local operators which can be written in the general form: 



where r nl and T n2 denote a combination of gamma matrices and q the quark flavor. They 
should respect the Dirac structure, the color structure and the types of quarks relevant 
for the decay being studied. They can be divided into two classes according to topology: 
tree operators (0±, O2), and penguin operators (O3 to 0\q). For tree contributions (where 
W ± is exchanged), the Feynman diagram is shown in Fig. |2] (left). The current-current 
operators related to the tree diagram are the following: 



where a and (3 are the color indices. The penguin terms can be divided into two sets. 
The first is from the QCD penguin diagrams where gluons are exchanged, while the 
second is from the electroweak penguin diagrams (where either a 7 or a Z° is exchanged). 
The Feynman diagram for the QCD penguin diagram is shown in Fig. |2] (right) and the 
corresponding operators are written as follows: 



Oi = (g^x g,,) (g fe r„ 2 gz) 



(30) 



0{ = qalti(l ~ 75)u0S/?7 M (l - 7b)6. 
O s 2 = qin{l - 7s)ms7^(1 - 75)6 , 



(31) 



3 



^(1-75)6^^(1-75)?', 




(32) 



5 



g 7M (l-7 5 )6^gV(l + 75)g' , 




(33) 



q' 



where q' = u, d, s, c. Finally, the electroweak penguin operators arise from the two Feyn- 
man diagrams represented in Fig. [3] (left) and Fig. [3] (right) where Z, 7 is exchanged from 



a quark line and from the W line, respectively. They have the following expressions: 



o 7 = 


2 TO 1 - 


-75)&^e^gV(l + 7 5 )g' , 




o 8 = 


3 


- 75)&/3 V^7 M (1 + 75)?a , 

q' 




o 9 = 


2^(1 - 


-7 5 )6^e g /gV(l-75)g' , 




o w = 


2^7m(! 


-75)6^^65*^(1-75)^ , 
<?' 


(34) 



where e q i denotes the electric charge of q'. 



4.1.2 Wilson coefficients 



As we mentioned in the preceding subsection, the Wilson coefficients [23J, Ci(/z), represent 
the physical contributions from scales higher than \i (the OPE describes physics for scales 
lower than \x). Since QCD has the property of asymptotic freedom, they can be calculated 
in perturbation theory. The Wilson coefficients include the contributions of all heavy 
particles, such as the top quark, the W bosons, and the charged Higgs boson. Usually, the 
scale \x is chosen to be of 0(m&) for B decays. The Wilson coefficients have been calculated 
to next-to-leading order (NLO). The evolution of C(/i) (the matrix that includes Cj(//)) 
is given by, 

C(ji) = U(n,M w )C{M w ) , (35) 
where U(/i, Mw) is the QCD evolution matrix: 



U(ji,M w ) 



An 



-J 



[/ u (/i,M, 



w , 



i a s {M w ) j 
Arc 



(36) 



with J the matrix summarizing the next-to- leading order corrections and U°(fi, Mw) the 
evolution matrix in the leading-logarithm approximation. Since the strong interaction 
is independent of quark flavor, the C(fi) are the same for all B decays. At the scale 
fj, = m b = 5 GeV, C(/i) take the values summarized in Table [T] To be consistent, the 
matrix elements of the operators, Oj, should also be renormalized to the one- loop order. 
This results in the effective Wilson coefficients, which satisfy the constraint, 



tree 



(37) 



here (Oi) tree are the matrix elements at the tree level. These matrix elements will be 
evaluated within the factorization approach. From Eq. (|3*Tj). the relations between C[ and 
Q are j21E5]: 



a 



Ci 
C 3 



Ps/S 

C 5 — C5 — P s /3 
Cj = C7 + P e , 



Cn — Ca 



Pe 



c 2 

C'e 



C 2 , 

C4 + P s , 



P., 



(38) 



where, 



P s = (« s /87r)C 2 (10/9 + G{m c , fi, q 2 )) 



P, 



and 



a em /97r)(3C 1 + C 2 )(10/9 + G(m e ,fi, q 2 )) 

m 2 — x(l — x)q 2 



G(m c , fx, q 2 ) = 4 / dx x(x — l)ln- 
'o 



(39) 



(40) 



Here q is the typical momentum transfer of the gluon or photon in the penguin diagrams 
and G(m c , fi,q 2 ) has the following explicit expression [26], 



&e G 



JmG 



2 
3 



.mi 5 A mz I ml 



1 - 4^1n 

« l-Jl 



4 n| 
9 2 




(41) 



Based on simple arguments at the quark level, the value of q 2 is chosen in the range 
0.3 < q 2 /m 2 < 0.5 OIHZl- From Eqs. (j38H41|) we can obtain numerical values for C[. 
These values are listed in Table 121 where we have taken a s {m z ) = 0.112, a em (m b ) = 
1/132.2, m b = 5 GeV, and m c = 1.35 GeV. 

4.1.3 Effective Hamiltonian 

In any phenomenological treatment of the weak decays of hadrons, the starting point is 
the weak effective Hamiltonian at low energy j2H]- It is obtained by integrating out the 
heavy fields (e.g. the top quark, W and Z bosons) from the standard model Lagrangian. 
It can be written as: 

^// = 7§E VcKAtCiMOiOi) , (42) 

* i 

where Gp is the Fermi constant, Vckm is the CKM matrix element (see Section I4~3|) . 
Ci(fjL) are the Wilson coefficients (see Section |4.1.2|) . Oj(/i) are the operators from the 
operator product expansion (see Section I4.1.1|) . and fi represents the renormalization 
scale. We emphasize that the amplitude corresponding to the effective Hamiltonian for 
a given decay is independent of the scale \i. In the present case, since we analyze direct 
CP violation in B decays, we take into account both tree and penguin diagrams. For the 
penguin diagrams, we include all operators O3 to Oio- Therefore, the effective Hamiltonian 
used will be, 



H eff = — 



G p 



10 



VuMdOl + C 2 O s 2 ) - V tb V* CiOi 



i=3 



H.c. 



(43) 



and consequently, the decay amplitude can be expressed as follows, 



A{B 



V ub V: s (C l {V l V 2 \O s 1 \B) + C 2 (V l V 2 \O s 2 \B))- 

10 

VnVZ^CiiVtVtlOilB) 



i=3 



+ H.c. , (44) 



where (ViV2\0i\B) are the hadronic matrix elements. They describe the transition be- 
tween the initial state and the final state for scales lower than /i and include, up to now, 
the main uncertainties in the calculation since they involve non-perturbative effects. 



4.2 New expression of helicity amplitudes hij according to Wil- 
son Coefficients 

4.2.1 General helicity amplitude 

The weak hadronic matrix element is expressed as the sum of three helicity matrix ele- 
ments, each of which takes the form H\(B — ► p°(o;)V2) = (V\V2\H w e ^ \B) , and is defined 
by gathering all the Wilson coefficients of both the tree and penguin operators. Linear 
combinations of those coefficients arise, such as (tree diagram contribution) and cj* 
(penguin diagram contribution). Then, in the case of B — > p°(u))V2, {V\ = p° or u), the 
helicity amplitude H\(B — > p°(uj)V2) has the general following expression: 

H X {B -> p\u)V 2 ) = (v ub V: s c p tl - V tb V t y pi ) Us aP , s e^(X)e*f(X)P^ 2 

+ i((3 p 2 e* V2 (X)e;(X) - (3 p 3 (e*y 2 (X).P B )(e;(X)-P B )) } + (k & K* s < - V tb V*c? 
Ple aM5 e; a {X)et{X)PlP 5 p + i(i3 p 5 e;(X)e* V2 (X) - P P (e;(X)-P B )(e*y 2 (X)-P B 

V ub V: s cl - V»V£<&) ^e a ^(X)e*f(X)P^P ] 



l>2 



s 

v 2 



s p - m l) + im u T t 

+ i(^e V2 (X)et(X) -^{e V2 {X)-P B ){el{X).P B )) j + [V ub V* s c^ 2 - V tb V t * s (% 2 

/3 4 %*/^C(A)4f(A)P^ , (45) 

with €v 2 , p ,oj(X) being the K*, p° and uj polarization vectors expressed in the B rest frame. 
Finally e a p^ is the antisymmetric tensor in Minkowski space. 

In Eq. ()45|) the parameters fa are mainly the form factors describing transitions be- 
tween vector mesons. They take the form: 

#i = ^fru*rn Vu *—£- V B ^{m\ y2 ) , (46) 

z m B + my 2 y 1 

G 

$i% = -&fv u vfnv u Vz{m B + m V2yi )A B ^ V2 ' Vl (m Vi y 2 ) , (47) 

= ^ffvuV 2 m Vl , V2 A B - V ^^ ) , (48) 

l m B + my 2 y 1 

here fv l ,v 2 ls either the p°,u or K* decay constant. V B ^ V2 ' Vl and A i ~ i 2 ' 1 are respec- 
tively the vector and axial form factors defined in Eqs. (|14H16j) . It is worth noticing that 
the tensorial terms which enter H\^B p Q (uj)V2) become simplified in the B rest frame 



because the four-momentum of the B is given by Pb = (tub, 0). Then, using the orthog- 
onality properties of e^(A), the helicity amplitude H\(B — > p°(a;)V2) acquires a much 
simpler expression than above: 



h x (b - P °(uj)v 2 ) = iB p x (v ub v: s( f tl - v tb v t y pi ) + ic p x {v ub v: s c p t2 - v tb v* 



P2> 



n 



P : ^J 



[s - m 2 ) + imj: u 



iB«{v ub v: s ci - v tb v* c « pi ) + icx(v ub v: s c? 2 - v tb v; s d : 



P2> 



(49) 



where the terms B^ 1 and C^ 1 take the following forms for the helicity (A) values, —1, 0, +1: 



B 



aVi m 2 B - {ml 2 + m 2 Vl ) _ „ \p\ 2 m 2 B 

A=0 - Pi r. ™ ^3 _. _. ■> 



Vi 



2 

riVi qV\ B 

°A=0 _ ^5 



2m V2 m Vl 

l V 2 + "M j 



m V2 m Vl 



K + ™v,) RVl \p\ 2 m 2 

h>6 



2m V2 m Vl " mv 2 mv 1 



(50) 

(51) 

(52) 
(53) 



In the above equations, |p| is the momentum common to the V\ and V2 particles in the B 
rest frame. It takes the form: 



\P\ 



\/[rn 2 B - {m Vl +my 2 ) 2 ][m| - (m Vl - m V2 ) 



2m 



(54) 



B 



where m\ and m 2 are the vector masses. Taking into account the previous relations, we 
arrive at the final form for the amplitudes H\\B — > p°(u)V2\: 

H x= o ±i (B - p°(u)V 2 



AX 2 



R 1 B._ + R 2 C._ 

A— ±1 * A— ±; 

n, 



rP dp 1 rP/^P 



+ 



s p - m 2 ) + imj, 



AA 



1 -qui /~iu) 

Xli-D\_0 T -rtobiO 
1 A— ±1 ^ A— ±1 



+ i 



ro; 1 TUJ/~1UJ 



where one defines, 



= rjX 2 ^ 1 - ym(c Vl 



Pi ) ' 



rVi 



pX 2 $+Me(c v p j) , 



(55) 



(56) 
(57) 



with V\ being either p° or uo. From Eq. (}55|) . the density-matrix elements /t A ,A' can be 
derived automatically and on has: 

h x>y = H X (B -> p\uj)V 2 )H* x ,{B -> p°MK 2 ) . (58) 

Because of the hermiticity of the matrix (/ia,A')> only six elements must be calculated. 



4.2.2 Explicit amplitudes for the B decays investigated 

By applying the formalism described in Section III, one gets in the case of the p° produc- 
tion, the following linear combinations of the effective Wilson coefficients: 
for the decay B° -> K*°p°: 



c 



' -n' + 9l r p - 3 (r 1 i C ' w i C i c * 

C'n \ 1 / „, . Cn 



4 = 0, < = _( Cl + _lj + _^ !o + _lj, (59) 

where are listed in Table 2. The coefficients, c£, relate to the tree diagrams and c£. to 
the penguin diagrams. To simplify the formulas we used N c for N^'* in the expressions 
(Eqs. (Hi-il). 
for the decay B~ ^ K*~ p°: 



hi 



< =c, * + w e > c -=^ + f + ^ + |- (60) 

In the case of u production one obtains the following linear combinations of effective 
Wilson coefficients: 
for the decay B° -> K*°u: 

0, cl=-C'-^ + Uc[ 1 ^ 



4 = c! + |, <ft = 2(c 5 + | + cj + |) + i(c? + f + C4 + |). (6i) 



for the decay i? — > cj: 

C{ C" / C" 

C + . = Cn + — — , C„, = C 4 + — — + Ci n + "77* 



We refer to Appendix A for details of the helicity amplitudes, while for the channel 
B ± — > p°(u)p ± we refer to Appendix B. 



4.3 CKM matrix and form factors 

In phenomenological applications, the widely used representation of the CKM matrix is 
the Wolfenstein parametrization m this approach, the four independent parameters 
are A c , A, p and 77. Then, by expanding each element of the matrix as a power series in 
the parameter A c = sin# c = 0.2209 (6 C is the Gell-Mann-Levy-Cabibbo angle), one finds 
(O(Xf) is neglected) 

/ 1-|A2 A c A\l(p- iV )\ 

Vckm = -A c 1 - \X 2 C AXl , (63) 

\A\Hl-p-ir,) -AXl 1 / 



where r\ plays the role of the CP- violating phase. In this parametrization, even though 
it is an approximation in A c , the CKM matrix satisfies unitarity exactly, which means, 



V CKM ' V CKM - I - VcKM ' VcKM • (64) 

The form factors, V(k 2 ) and Aj(k 2 ), depend on the inner structure of the hadrons. Here we 
will adopt two different theoretical approaches. The first was proposed by Bauer, Stech, 
and Wirbel ^H] (BSW), who used the overlap integrals of wave functions in order to 
evaluate the meson-meson matrix elements of the corresponding current. In that case the 
momentum dependence of the form factors is based on a single-pole ansatz. The second 
approach was developed by Guo and Huang (GH) [29J, who modified the BSW model by 
using some wave functions described in the light-cone framework. Nevertheless, both of 
these models use phenomeno logical form factors which are parametrized by making the 
assumption of nearest pole dominance. The explicit k 2 dependence of the form factor 

is mum]: 



v(k 2 ) = -t^-tt , Mk 2 ) = . N , (65) 



k' 2 



k 2 

5] 



where . and my are the pole masses associated with the transition current and hy and 
Ha. are the values of the form factors at q 2 = 0. 

5 Monte-Carlo simulations: computation of hij and 
general results 

5.1 Numerical inputs 

5.1.1 CKM values 

In our numerical calculations we have several parameters: q 2 , N%** and the CKM matrix 
elements in the Wolfenstein parametrization. As mentioned in Section IV, the value of q 2 
is conventionally chosen to be in the range 0.3 < q 2 /m b 2 < 0.5. The CKM matrix, which 
should be determined from experimental data, is expressed in terms of the Wolfenstein 
parameters, A, A c , p, and rj |T9j . Here, we shall use the latest values [3*T] . which were 
extracted from charmless semileptonic B decays, (|Vu&|), charmed semileptonic B decays, 
(|V^,|), s and d mass oscillations, Am s , Am,i, and CP violation in the kaon system (e^), 
(p,T]). Hence, one has, 

A c = 0.2237 , A = 0.8113 , 0.190 < p < 0.268 , 0.284 < rj < 0.366 . (66) 

These values respect the unitarity triangle as well (see also Table |SJ). In our numerical 
simulations, we will use the average values of p and r\. 

5.1.2 Quark masses 

The running quark masses are used in order to calculate the matrix elements of penguin 
operators. The quark mass is evaluated at the scale p ~ m;, in B decays. Therefore one 



has HH, 



m u (p = m h ) = 2.3 MeV , fn d (p = m b ) = 4.6 MeV , 

m s (n = m h ) = 90 MeV , m b (/i = m h ) = 4.9 GeV , (67) 

which corresponds to m s (fi = 1 GeV) = 140 MeV. For meson masses, we shall use the 
following values [T3] : 

m B ± = 5.279 GeV , m K * = 0.896 GeV , m u = 0.782 GeV , 

m B o = 5.279 GeV , m p ± = 0.770 GeV , m^± = 0.139 GeV , 

mjfri = 0.892 GeV , m p o = 0.770 GeV , m^o = 0.135 GeV . (68) 



5.1.3 Form factors and decay constants 

In Table E] we list the relevant form factor values at zero momentum transfer [TTH EHJ EHJ 
for the B — > if*, B —>■ p and B uj transitions. The different models are defined 
as follows : model (1) is the BSW model where the q 2 dependence of the form factors 
is described by a single- ansatz. Model (2) is the GH model with the same momentum 
dependence as model (1). Finally, we define the decay constant for vector (_/y) meson as 
usual by, 

qil^lfy = f P m p ep for p and otherwise , 
(V(q)\qa,q 2 \0) = f v m V ev , (69) 

with q being the momentum of the vector meson and my and ey being the mass and 
polarization vector of the vector meson, respectively. Numerically, in our calculations, we 
take [13], 

f K * = 214 MeV , f p = 221 MeV , f u = 195 MeV . (70) 

Finally, the free parameter (effective number of color, N^-') is taken to lie between the 
lower(upper) limits 0.66(2.84) for b — > s transition. Nevertheless, we focus our analysis 
on values of N£** bigger than 1, as suggested in [33]. Regarding the b — > d transition, the 
lower(upper) limits for N e Jf are 0.98(2.01) [H]. 



5.2 Simulation of the p° — uj mixing 

All the channels studied here include at least one p° meson which mixes with the u meson. 
The other vector mesons are either a K* 0,± or a p^. Thus, the mass of each resonance is 
generated according to a relativistic Breit-Wigner: 



da TrM r 
dM 2 N (Mi-Mlf + (T R M R 



^777^ 77772 , 77 77 x 2 ' ( 71 ) 



where is a normalization constant. In Eq. (|71|). Mr and T R are respectively the mass 
and the width of the vector meson which have been determined experimentally. M is 
the mass of the generated resonance. A simple and phenomenological relation describing 



the amplitude for p° — u mixing is used for the Monte-Carlo simulations 35 j. In the 
expression for the Breit-Wigner, the p°-propagator is replaced by the following one: 

1 1 T n 

— = - + ^^, (72) 

where and T p are respectively the uj and p production amplitudes. In addition, n pw 
is the mixing parameter for which recent values come from e + e _ — > tt + tt~ annihilations. 
Explanations have been already been given in Section III. Finally, 1/sy has the same 
definition as in Eq. (|2*2*j). Because the same physical processes enter the production of 
both the p° and u resonances (they are both made out from uu and dd quark pairs with 
the same weight 1/2), it seems natural to choose T w /T p = 1. So, the invariant mass 
distribution of the 7r + 7r~ system becomes simplified, being given by, 

da/dm 2 oc\A(p\uj))\ 2 , (73) 

where A(p°(u>)) is the amplitude of the two mixed Breit-Wigner distributions and m is 
the 7r + 7r~ invariant mass. In Fig. IH the 7r + 7r~ invariant mass spectra for p° — u mixing 
is displayed. Because of the very narrow width of the u, (r w = 8.44 MeV), we notice a 
high and narrow peak at the u pole (~ 782 MeV). 



5.3 Density matrix h xx > 

Three main parameters remain free in our simulations: the ratio q 2 /m 2 (related to the 
penguin diagrams), the form factor model (GH or BSW) and the effective number of 
colors, N^** (associated with the factorization hypothesis). The histograms plotted in 
Fig. El display spectra of the diagonal and normalized density matrix elements hi^, for 
the channels B° — > p°(uj)K*° (left hand-side) and B + — > p°(uj)p + (right hand-side). The 
input numerical parameters are q 2 /m 2 b = 0.3, N%** = 2.84 (left hand-side figure) or 
jyeff = 2.01 (right hand-side figure), and the GH form factor model is applied for both 
decays. Note also that the average values of CKM parameters p and r\ are used. The wide 
spectrum of values of the density matrix element h\ } \, is caused by the resonance widths 
(especially that of the p) which provides, in turn, a large spectrum for the common 
momentum py in the B rest-frame. Whatever the p (u)V2 channel is, h 00 = |-f^o| 2 ) 
which corresponds to longitudinal polarization, is the dominant value. Numerically, for 
the B° — > p°(u)K*° decay, the mean value of h 00 is around 0.87 while it is of order 
0.90 for B + — > p°(uj)p + . The dominance of the longitudinal polarization has already 
been confirmed experimentally, since recent experimental data related to the channel 
B — » J / ipK* show clearly that the longitudinal decay amplitude dominates in that case, 
with \Hq^ = 0.59 ± 0.06 ± 0.01 Extrapolating these results to the charmless vector 
meson final states requires some modifications of the form factors without a big change of 

I 1 2 

the relative contributions of the polarization states. Regarding h = , it represents 

less than 0.5% of the total amplitude for both decays. This numerical result is confirmed 
by complete analytical calculations. 

In Figs. Eland the real and imaginary parts of the non-diagonal and normalized den- 
sity matrix elements hij are shown for the channels B° — > p (u)K* and B + — > p°(o;)p + , 
respectively. The input parameters are the same as previously mentioned. The main 
feature of the non-diagonal matrix elements, hi^, is the smallness of both the imaginary 



and real parts - the imaginary part being at least one order of magnitude smaller than 
the real part one. For the B + — ► p°(u))p + decay, we observe that the mean value of all 
the imaginary parts is zero, whereas it can vary for the other decay. Note also that each 
of the three real parts are quite similar for both decays. Because of the tiny value of 

h = |-ff_i| , the moduli of the non-diagonal elements, h + _ = H + H^_ and h _ = H H^, 

are very small, while the modulus of h +0 = H + Hq is around 0.3 for both decays. As a 
first conclusion, the general behavior of the density matrix seems to be similar whatever 
the decay is. Experimentally, only the mean values of the diagonal elements and will 
be able to be measured through the angular distributions. 

These angular distributions are plotted in Figs.|S]and|n]in the helicity frame and in the 
transversity frame, respectively for B° — > p°(u)K*° and for the usual input parameters. 
Their normalized pdfs have been displayed above in Eq. (JTUJ). As a consequence of the 
small value of (h + _), the azimuthal angle distribution in the helicity frame is nearly flat, 
whereas it is sinusoidal in the transversity frame. From the distribution as a function 
of polar angle (in the TF) displayed in Eq. (|13p. one can infer a mean value of the 
amplitude. This represents an additional piece of information through which one can 
access the dynamics of B(B) decays into two charmless hadrons. 

6 Branching ratio and asymmetry in B decays into 
two vector mesons 

The analytic expressions for the density matrix elements, allow us to calculate the 
hadronic branching ratios SSiB — > p°{uj)V2) and to estimate the asymmetries related to 
B and B decays. All these physical observables depend primarily on a subset of the 
parameters mentioned previously, such as the form factors, the ratio q 2 jm 2 (where q 2 
is the mass of the virtual gluon in the penguin diagram), the effective number of colors, 
Nl*' (used as a free parameter in the framework of the factorization hypothesis), and the 
CKM matrix element parameters p and r\. 

6.1 Branching ratio: results and discussions 

Departing from the definition of the branching ratio {3B{B —►/)), 

*< B ^> = TTrai)< (74 > 

the width F(B — ► V1V2) can be inferred from its differential form given by the standard 
relation [3*7j : 

<^( B -™) = ^^(B^v^^-^-sHP-n-n). (75) 

In Eq. (|75|) . P = (M, 0), where M = and (E\,pi) and (E^jpz) are the 4-momenta 
of V\ and V2, respectively, in the B rest frame. Because of the large width of the p° 
meson (T p m 150 MeV) and the K* meson (r^» w 50 MeV), the energy, Ei, and the 
momentum, pi, of each vector meson vary according to the generated event. Computa- 
tion of T(B — > p°(a>)V2) could not be done analytically but numerically by Monte-Carlo 



methods. A total number of 50000 events have been generated in order to obtain a precise 
estimate of this decay width. 

In Tables El and El we list (respectively) the branching ratios for B — > p°(u)V2 and 
B — ► p°{uj)V2 and their dependence on the form factor models (BSW and GH), q 2 /ml, 
iV c e// and the average values of the CKM parameters p and rj. For a fixed value of q 2 jm 2 , 
there are important variations of the branching ratios, depending on the form-factor 
model. They can vary by up to a factor two. In the framework of a given form-factor 
model, some branching ratio modifications appear with q 2 /rrib 2 , especially in the channels 
including a K* . However, these changes do not exceed 34%. Regarding the ratio between 
&(B° -> p°(uj)K*°) and M(B+ -> P °(uj)K* + ), its value is found to be of the order 0.40 
for the BSW model and 0.34 for the GH model. 

Finally, we observe that the relative difference between two conjugate branching ratios, 
SS{B — > /) and SS[B — > /), is almost independent of the form-factor models, for a fixed 
value of q 2 /m b 2 . It can be computed from the two tables just mentioned and, usually, it 
does not exceed 20%. The exception is for the i^* ± p°(u;) channels, where it reaches 39%. 



6.2 Asymmetry: results and discussions 

A search for direct CP violation requires asymmetries between conjugate final states 
coming from B and B decays respectively. In our case, these searches are performed 
in two complementary ways. We consider first the global CP-violating asymmetry acp, 
calculated from branching ratios: 

&{B -> f) - WjB -> f) 

acp = _ — . (7o) 

9{B - /) + S9(B - /) 

Secondly, we use the partial widths of B(— > /) and B{— > /), calculated as described 
above together with the differential asymmetries investigated as a function of the 7r + 7r~ 
invariant mass in the whole range of the p° Breit-Wigner resonance. Hence, acp(m) takes 
the following form: 

r m (B ->/) - rug -> /) 

aCF(m) = T m (B^f) + T m (B^f) • < 77 > 

where m is the ir + n~ invariant mass. T m (B f) and T m (B —* /) in Eq. (|77j) are the 
partial widths written as a function of m. 

In Tabled we list the global CP-violating asymmetry between the B and B decays 
for the channels under investigation. It can be noticed that, for a fixed value of q 2 Im 2 , 
the two form factor models provide quite similar results. For different q 2 jm 2 values, the 
corresponding results could vary, especially in the i\"* ± p°(o;) channels. In Figs. 1101 and fTTI 
we show, respectively, the histogram of the direct CP- violating asymmetry parameter 
a cp{ m )i f° r the decays B° — > p°(u)K*° and B + — > p (u)p + , as a function of the tt + tt~ 
invariant mass in the u mass region and for both form factor models. The asymmetry 
reaches its maximum when y/s is around 780 MeV. However, outside the displayed win- 
dows, the asymmetry goes to zero in any case. The peak of the asymmetry is emphasized 
when the GH form factor model is used in our simulations. For the K*°p°(uj) channels, 
the maximum of the CP violating asymmetry is around 13% and 16%, for the BSW model 
and the GH model, respectively. Finally, we emphasise that the p ± p°(uj) channels present 



the most intriguing results because, in any case, their asymmetry is at least 80% (BSW 
model) and can reach 95% (GH model). This last channel is highly recommended for a 
direct search for CP violation. 

7 Perspectives and conclusions 

We have studied direct CP violation in decay process such as B — > p°(o>)V2 — > 7i + n^V2, 
where Vi is either K* 0,± or p^, with the inclusion of p° — uj mixing. When the invariant 
mass of the 7t + tt~ pair is in the vicinity of the u resonance, it is found that the CP- 
violating asymmetry, acp(m), reaches its maximum value. In our analysis we have also 
investigated the branching ratios for the same channels. Thanks to the standard helicity 
and transversity formalisms, rigorous and detailed calculations of the B 0± decays into two 
charmless vector mesons have been carried out completely. Using the effective Hamilto- 
nian based on the operator product expansion with the appropriate Wilson coefficients, 
we derived in detail the amplitudes corresponding to B — > p (u)V2 — > n + n^V2 decay and 
the density matrix, h\y as well. 

In order to apply our formalism, we used a Monte-Carlo method for all the numer- 
ical simulations. Moreover, we dealt at length with the uncertainties coming from the 
input parameters. In particular, these include the Cabibbo-Kobayashi-Maskawa matrix 
element parameters, p and r), the effective number of colors, A^f , coming from the naive 
factorization and two phenomenological models in order to show the possible dependence 
on form factors, GH or BSW. These form factors vary slightly according to the final 
states. Recall that this work was achieved by applying a phenomenological treatment, 
where some assumptions regarding the evaluation of the hadronic matrix elements have 
been made. In this approach, corrections associated with the limit of validity of the 
factorization hypothesis were parameterized phenomenologically and may involve large 
uncertainties. 

As a major result, the predominance of the longitudinal polarization, /i o, has been 
pointed out in all the investigated decays. We also found a large direct CP- violating 
asymmetry in these B decays into two charmless vector mesons. We stress that, without 
the inclusion of p°—u) mixing, we would not have a large CP-violating asymmetry. Finally, 
we predicted branching ratios to be of the order 0.7 — 2.1 x 10~ 6 for K*°p°(u) and of the 
order 2.3 — 6.6 x 10 -6 for A'* ± p°(o;) (depending on the different phenomenological models). 
For the channel p ± p°(oj), we found the branching ratios to be of the order 11 — 24 x 10~ 6 . 

Two main conclusions can be drawn. The first is the relative importance of the form 
factor model which is used, since some branching ratios in B — > p°(u))V2 could change by 
up to a factor two. The second is the important role of p° — uj mixing, which can enhance 
considerably the asymmetry parameter acp, between the conjugate final states coming, 
respectively, from B and B decays. 

Beside the "standard" ways to look for direct CP violation, such as the difference 
between branching ratios and/or discrepancies in the angular distributions of the decay 
products, we have presented a detailed discussion of a new method. This involves the 
variation of acp as a function of the invariant mass over the whole range of the p° 

resonance [T2*l l34|. We believe that this method will be very fruitful for future experiments 
and has already been implemented in the generator of the LHCb experiment. Indeed, we 
look forward to being able to apply the formalism developed here to the analysis of 



experimental data for decays such as B — > p°(o;)X^ (with V2 being either a if* or a p ± ) 
in the near future. 
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Appendix 

A Practical calculations of the helicity amplitudes 

The helicity formalism in the case of vector mesons requires the introduction of three 
polarization four-vectors for each spin 1 particle 



e(l) = (0,e(l)) , e(2) = (0,e(2)) , and e(3) = (\k\/m,Ek/m) . (78) 

They also satisfy the following relations as well: 

e(if = -1 , and e(i) • e(j) = , with i^j, (79) 

where m, E and k are respectively the mass, the energy and the momentum of the vector 
meson, k is defined as the unit vector along the vector momentum, k = k/\k\. The 
three vectors e(l),e\2) and <f(3) = Ek/m form an orthogonal basis. e(l) and e*(2) are the 
transverse polarization vectors while e(3) is the longitudinal polarization vector. These 
three vectors allow one to define the helicity basis: 

(e(l) + u(2)) (e(l)-ie(2)) 
e(+) = -j= , e(-) = -j= , and e(0) = e(3) . (80) 

These 4-vectors are eigenvectors of the helicity operator TC corresponding, respectively, 
to the eigenvalues A = +1,-1 and 0. In the B 0± rest-frame, the vector mesons have 
opposite momentum k\ = —k 2 and their respective polarization vectors are correlated. 
This implies the following expressions, 

(k sin 9 cos 1 
k sin 9 sin < 
k cos 9 

where 9 and are respectively polar and azimuthal angles of the produced K*. In our 
case, one has for the transversal polarization vectors (K* and p) the expressions: 

cos 9 cos 
ek (1) = I cos#sin0 I — e p (l) 
— sin 9 



and, 



- sm q 

e K {2) = | cos0 | = -e p (2) 

Regarding the longitudinal polarization, e^(3) and e p (3) take the form: 

£jf (3) = [ A, ^ , e p (3) = f ^(-^ . (81) 
m K m K J \m p m p J 

By applying the relations from Eq. (|8U|) . one can expressed vectors e{i) in the helicity 
basis and one gets e*(±): 



I cos 9 cos — i sin i 
cos 9 sin + i cos ( 
y — sin 9 



(82) 



= e p (+) . (83) 

The weak hadronic amplitude is therefore decomposed, in the helicity basis, according to 
the general method developed by Bauer, Stech and Wirbel JT0| - This will allow one to 
obtain two interesting results. Firstly, one can isolate the contribution of each helicity 
state to the total amplitude. Secondly, the contributions of the tree and penguin operators 
to the total amplitude can be separated via the helicity states. 

The knowledge of the main input parameters p, r], A, sin# c (= A c ) and the masses and 
widths of the intermediate resonances allow a complete determination of the three helicity 
amplitudes H\(B — > p G {oj)V<2), where the helicity A can take the values -1,0 or +1. 



cos 9 cos + i sin 
cos 9 sin — % cos 
— sin 9 



/V2 



B Channel B ± p°(oj)p ± 

The formalism applied in case of B — > p°(uj)K* can be extend to B ± — >• p (uj)p ± . Nev- 
ertheless, in the last case one has b — > d transition instead of b — > s. The amplitude 
H\(B — > p°{uS)V2) has the form: 



H x=0±i (B - p°(u;)V 2 ) 



AX 3 



RiB,_ + RnC o 
A— ±1 ^ A— ±1 



+ i 



i l- D A = q . + J 2°A=° - 



n 



(s p - ml) + im u Y u 



AX' 



+ Z 



1 \ a \=° + h °A=0- 



54) 



where one defines, 



2 



= (1 - + v - (1 - p) ym(c^) , (85) 



2 



= (1 - y )pcf + V Sm(<%) + (1 " p) , (86) 

with V\ being either p° orw. 

If V 1 = p and i = 2 then i^ 1 = /f 1 = . (87) 

The expressions for c^ 1 and c^ 1 , which correspond to the investigated channel, take the 

following form: 

for the decay B~ — > p°p~ : 



41 iV c N c 



In the case of oj production, one obtains the linear combinations of the effective Wilson 
coefficients: 

for the decay B~ — > up~: 



c 



C 1 C / C 1 

ii I ' 1 Mi _ r<i i 3 i u 9 



" ~ c ' 2 + jv; • c '>- c ' + it + K ' 1 " + jv^J • 

+ ^ + f + ^ + f- C ; -g). (89) 

All the terms used in the appendix have been defined in Section IV. 
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Fig. E Transversity frame for B — > p°K*. 

Fig. 121 Tree diagram (left), and QCD-penguin diagram (right), for B decays. 

Fig. El Electroweak-penguin diagram (left), and electroweak-penguin diagram with 
coupling between Z, 7 and W (right), for B decays. 

Fig. 0] Spectrum of p° — uj mixing (in MeV/ c 2 ) , simulated by the interference of two 
Breit-Wigner curves. 

Fig. Spectrum of h , h 00 , h ++ . Histograms on the left correspond to the channel 

B° — > p°(lv)K*° where the parameters used are: q 2 /m 2 b = 0.3, N^* = 2.84, p = 
0.229, 1] = 0.325 and form factors from the GH model . Histograms on the right 
correspond to the channel B + — > p°(uj)p + for the same parameters with N%" = 2.01. 

Fig. El Spectrum of ^eih^) and J ! m{hij) where % 7^ j. Histograms correspond to the 
channel B° — > p°(u)K*° where the parameters used are: q 2 /m 2 b = 0.3, N%** = 2.84, 
p = 0.229, 7] = 0.325 and form factors from the GH model. 

Fig. [3 Spectrum of £%k(hij) and J ! m{hij) where i 7^ j. Histograms correspond to the 
channel B + — > p°(u)p + where the used parameters are: q 2 /m 2 b = 0.3, N%" = 2.01, 
p = 0.229, 7] = 0.325 and form factors from the GH model. 

Fig. |S1 Spectrum of polar angle (upper figure) and azimuthal angle (lower one) in 
the helicity frame for the channel B° p (u)K*°. Parameters are: q 2 /m 2 b = 0.3, 
iV c e// = 2.84, p = 0.229, r] = 0.325 and form factors from the GH model. 

Fig. |H1 Spectrum of polar angle (upper figure) and azimuthal angle (lower one) in the 
transversity frame for the channel B° — > p°(uj)K*° . Parameters are: q 2 /m 2 b = 0.3, 
iV c e// = 2.84, p = 0.229, r] = 0.325 and form factors from the GH model. 

Fig. El CP-violating asymmetry parameter acp{m), as a function of the tc + 7t~ 
invariant mass in the vicinity of the uj mass region for the channel B — > 
Parameters are: q 2 /m 2 b = 0.3, iV c e// = 2.84, p = 0.229, rj = 0.325. Solid triangles 
up and circles correspond to the BSW and GH form factor models respectively. 

Fig. ^2 CP-violating asymmetry parameter acp{m), as a function of the tc + 7t~ 
invariant mass in the vicinity of the u mass region for the channel B + — > p°(uj)p + . 
Parameters are: q 2 /m 2 b = 0.3, N e J f = 2.01, p = 0.229, 77 = 0.325. Solid triangles 
down and circles correspond to the BSW and GH form factor models respectively. 
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• Tablets , B~ branching ratios (in units of 10~ 6 ) using either the BSW or GH form 
factor models, for q 2 /m 2 b = 0.3(0.5), with N^ x = 2.84(2.82), N b c ^ x = 2.01(1.95), 
p = 0.229 and 77 = 0.325. 

• TablelH-B , B + branching ratios (in units of 10~ 6 ) using either the BSW or GH form 
factor models, for q 2 /m 2 = 0.3(0.5), with N^ x = 2.84(2.82), N b c ^ x = 2.01(1.95), 
p = 0.229 and rj = 0.325. 

• Table [7| Global CP-violating asymmetries (in percents) using either the BSW or 
GH form factor models, for q 2 /m 2 b = 0.3(0.5), with N^ x = 2.84(2.82), N^ x = 
2.01(1.95), p = 0.229 and rj = 0.325. 




Figure 2: Tree diagram (left), and QCD-penguin diagram (right), for B decays. 



Figure 3: Electroweak-penguin diagram (left), and electroweak-penguin diagram with 
coupling between Z, 7 and W (right), for B decays. 




Figure 4: Spectrum of p° — to mixing (in MeV/c 2 ), simulated by the interference of two 
Breit-Wigner curves. 





Figure 5: Spectrum of h ,hoo,h ++ . Histograms on the left correspond to the chan- 
nel B° — > p°{ui)K*° where the parameters used are: q 2 /m 2 t > = 0.3, N^* = 2.84, 
p = 0.229,r/ = 0.325 and form factors from the GH model. Histograms on the right 
correspond to the channel B + — > p°(uo)p + for the same parameters with N^f = 2.01. 
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Figure 6: Spectrum of J^e(hij) and J^m(hij) where i ^ j. Histograms correspond to 
channel B° — > p°(u)K* where the used parameters are: g 2 /m 2 5 = 0.3, N^f = 2.84, 
p = 0.229, f] = 0.325 and form factors from the GH model. 




Figure 7: Spectrum of Me(hij) and J^m(hij) where % ^ j. Histograms correspond to 
the channel B + — > p°(uj)p + where the used parameters are: q 2 /m 2 ^ = 0.3, N^f = 2.01, 
p = 0.229, i] = 0.325 and form factors from the GH model. 
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Figure 8: Spectrum of polar angle (upper figure) and azimuthal angle (lower one) in the 
helicity frame for the channel B° — > p°(uo)K*°. Parameters used are: q 2 /m 2 b = 0.3, 
iV c e// = 2.84, p = 0.229, rj = 0.325 and form factors from the GH model. 




Figure 9: Spectrum of polar angle (upper figure) and azimuthal angle (lower one) in 
the transversity frame for the channel B° — > p°(ui)K*°. Parameters are: q 2 /m 2 b = 0.3, 
iV c e// = 2.84, p = 0.229, rj = 0.325 and form factors from the GH model. 
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Figure 10: CP-violating asymmetry parameter acp{m), as a function of the TT + ir invari- 
ant mass in the vicinity of the uj mass region for the channel B° — > p°(ui)K*°. Parameters 
are: q 2 /m 2 b = 0.3, iV c e// = 2.84, p = 0.229, rj = 0.325. Solid triangles up and circles 
correspond to the BSW and GH form factor models respectively. 
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Figure 11: CP- violating asymmetry parameter acp{m), as a function of the 7r + 7r invari- 
ant mass in the vicinity of the u mass region for the channel B + — > p°(uj)p + . Parameters 
are: q 2 /m 2 b = 0.3, iV c e// = 2.01, p = 0.229,7? = 0.325. Solid triangles down and circles 
correspond to the BSW and GH form factor models respectively. 



Ci(/i) for jj, = 5 GeV 



Ci 


-0.3125 






+ 1.1502 




C 3 +0.0174 




+0.0104 


C 4 -0.0373 


c 6 


-0.0459 


C 7 -1.050 x 10" 5 


c 9 


-0.0101 


C 8 +3.839 x 10" 4 


Cio 


+ 1.959 x 10~ 3 



Table 1: Wilson coefficients to the next-leading order (see the reference in text). 



c[ 


q 2 /ml = 0.3 




q*/ml = 0.5 






-0.3125 




-0.3125 




c> 


+ 1.1502 




+ 1.1502 






+2.433 x 10~ 2 + 1.543 x 10" 


-H 


+2.120 x 10~ 2 + 2.174 x 10" 


-H 




-5.808 x 10- 2 -4.628 x 10" 


-H 


-4.869 x 10~ 2 - 1.552 x 10" 


-H 




+1.733 x 10~ 2 + 1.543 x 10" 


-H 


+1.420 x 10~ 2 + 5.174 x 10" 


-H 




-6.668 x 10" 2 - 4.628 x 10" 


-H 


-5.729 x 10~ 2 - 1.552 x 10" 


-H 


c 7 


-1.435 x 10~ 4 - 2.963 x 10" 


-H 


-8.340 x 10~ 5 - 9.938 x 10" 


-H 




+3.839 x 10" 4 




+3.839 x 10" 4 






-1.023 x 10" 2 - 2.963 x 10" 


-H 


-1.017 x 10~ 2 - 9.938 x 10" 


-H 


C\Q 


+ 1.959 x 10~ 3 




+ 1.959 x 10" 3 





Table 2: Effective Wilson coefficients related to the tree operators, electroweak and QCD 
penguin operators (see the reference in text). 





a 


P 


7 


\Pmini ^]rain) 


104°47 


19°32 


56°21 


(Pmin i 1]max ) 


93°13 


24°31 


62°56 


(.Pmax i Vmin) 


112°14 


21°20 


46°66 


(Pmaxi Tjmax) 


99°66 


26°56 


53°78 



Table 3: Values of the CKM unitarity triangle for limiting values of the CKM matrix 
elements. 



B^p 







hy 


hA — ^Ai 


h Al 


h A2 


m v (GeV 2 ) 


m Ai (GeV 2 ) 


model 


(1) 


0.329 


0.281 


0.283 


0.283 


5.32 


5.32 


model 


(2) 


0.394 


0.345 


0.345 


0.345 


5.32 


5.32 


B -> LJ 






hy 


hA = h A 3 


h Al 


hA 2 


my (GeV 2 ) 


m Al (GeV 2 ) 


model 


(1) 


0.328 


0.280 


0.281 


0.281 


5.32 


5.32 


model 


(2) 


0.394 


0.345 


0.345 


0.345 


5.32 


5.32 


B ^ K* 






hy 


hA - h A 3 


h Al 


h A 2 


my (GeV 2 ) 


m Al (GeV 2 ) 


model 


(1) 


0.369 


0.321 


0.328 


0.331 


5.43 


5.43 


model 


(2) 


0.443 


0.360 


0.402 


0.416 


5.43 


5.43 



Table 4: Form factor values for B — > p, £> — > and i? — > X* at g 2 = (see the reference 
in text). 
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9 2 
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0.3 


2.1 


1.0 


B° -> K*°p°(uj) 










0.5 


1.5 


0.73 




0.3 


6.6 


3.9 


B~ -> K*-p°(u) 
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6.2 


3.6 
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24 
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Table 5: B°,B branching ratios (in units of 10 6 ) using either the BSW or GH form 
factor models, for g 2 /m 2 = 0.3(0.5), with N b c ^ x = 2.84(2.82), N b c ^ x = 2.01(1.95), p = 
0.229 and rj = 0.325. 
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GH 










0.3 


2.1 
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B° -> K*°p°(u) 
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1.7 


0.88 




0.3 


5.8 


3.4 


B+ -> K*+p (u) 










0.5 


3.8 


2.3 




0.3 


20 


11 












0.5 


20 


11 



Table 6: B°,B + branching ratios (in units of 10 6 ) using either the BSW or GH form 
factor models, for q 2 /m 2 b = 0.3(0.5), with N^ x = 2.84(2.82), N^ x = 2.01(1.95), p = 
0.229 and rj = 0.325. 





channel 


i 2 
™% 


BSW 


GH 


K* (K*°)p°(uj) 


0.3 
0.5 


+0.36 
+4.70 


-0.45 
+5.90 


K*- 




0.3 
0.5 


-6.6 
-23.0 


-6.37 
-22.0 




-(p + )pV) 


0.3 
0.5 


-8.5 
-8.7 


-9.6 
-9.9 



Table 7: Global CP-violating asymmetries (in percents) using either the BSW or GH 
form factor models, for q 2 /m 2 b = 0.3(0.5), with N^ x = 2.84(2.82), N^ x = 2.01(1.95), 
p = 0.229 and rj = 0.325. 



